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eapite treaendous technolojical discoveries and developmenta, 
the gas and ofl industry of teday ia still plagued with unresolved and 
unexplained problems ami phenomena, ospecially in the field of reservoir 
engineering. One of these ie the macroscopic flow of fluids in reservoir 
reck formations «here <ae and oil are found end produced, This thesis is 
@ stady of but a smail facet of the preblem and has been conducted with 
the idea that any work done in this ares will help contribate to an in- 
creased knowledge of the subject. 

The author wishes to seknowleige the sincere and helpfal suggea- 
tions, advice, and encouragenent of Professor 4olbrook 5. Deteet, Head of 
the Petrolews Engineering Sepertaent, University of Pittsburgh, in the 
Projection, progression, and successful eompletion of this study. 
Acknowledgment is also uade te those authors, investigators, and researchers 
whose works ware utilized for beckground and basis of the study without 
which such en undertaking would have been impoasible. 
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I. LETRODUCTION 


the Petroleum Industry finds and obtains ite raw materials, ces 
and ofl, fay beneath the surface of the earth in srese of rock formation 
called reservoirs, These reservoire ere not, as is usually presumed by 
the uninformed, large voids filled with gas and of] which, when tapped, 
say be easily produced. (m the contrary, they are formations of porous 
and permeable reek with email, uacroscepic, open spaces between the 
grains of the satrix saterial in shieh the gas and ofl is stored; a 
fagter which makes the production of such stored gas end oil difficult 
ami, in some cages, igpossible. it ie for this reasen that the complete 
stady of all conmiitions and probabilities of fluid flew in porous media 
becomes an iaportant facet in the preduction of gas and oil. 

The purpose of this stady war to investigate and attenpt to 
evaluate the inter-relationships of some of the major physical properties 
of porous reserveir rock which affeot the flow of fluids through it, 
Although those intererelationships have been expressed by several 4i f+ 
ferent formulations, mentioned later, ell designed to elucidete fluid 
flow, it was thought thet eo more usefal and informative modi fication 
could be developed and applied te the problen of correlating permeability, 
porosity and grain sise of porous aedia te fluid flow throagh such media, 
Sueh a study, designed te reveal the sechenics of aacroscopie fluid flew 
in porous media, should be usefal in evaluating some of the many factors 
whieh suet ultimately be known if 4 satisfactory solution of the probles 
ie ever to be obtained. Any adiitional knowledse of these factors should 
contribute to a more complete understanding of the complexities of gas 
and oil production from the reservoir, 
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Yo the petroleum engineer, especially one interested in reservoir 
enginoering, @ very important physical property of reservoir rock ie 
permeability, which can be properly defined ag a measure of the flaid- 
traneal iting capacity of & porque material, or as the abliity of ea fluid 
to flow within the interconnected pere spaces of & porous material. A 
second physical property, aleo equally important, is porveity. In oil] and 
g@e8 reservoirs, porosity represents the percentage ef total rock volume 
which is available for ocoupancy by either liquids or gueses or both. In 
connection with the tera porosity, care should be taken to distinguish 
veteeen absolute and effective porosity, Abselute poresity is the per- 
centage of void space with respect to total volume. Effective porosity, 
on the other hand, is the pereentage of intereonnected void space in the 
total volume, A third physical property of reservoir rock is actual 
grein sise andi grein sise distribution of the matrix saterial which fons 
the reservoir rock. 

Lt has been demonstrated by experiment and empirical formulas 
that permeability, effective poresity, and grain size are interrelated, 
ene depending upon the other to varying degrees. In porous material in 
which the pores are inter-connected, there exists no definite relation 
between permeability and effective porceity. ‘These, however, are directly 
influenced by grain size, in that grain size or particle size and ite 
attendant factors of uniformity of grain sise, shape of the grain, sanner 
of grein packing, and the amount of cementing material between the grains 
(degree of lithification) determine the sise and shape of the pore open- 
ings, the extent of interconnections of the pore spaces, and the ratio 
of their volume to the total volume of the rock. 
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have indleated that actual «rain sige 
and effective porosity are excellent criteria for judging permeabilities 
of various porous acdia, In mach of the literature, variations of 
perneability and porosity within individual reservoire have been found te 
be inter-related to a degree and have been shom by Pettke4s? in the case 
of the Bradford Third sand of the Sradford field of Pemeylvania to have 
these relationebipa 

K @ beu9 x 107% g 5065 
and Kw Sed x 1 g )**? 
where K is the permeability and @ is the porosity, The first forms was 
deterained from oarlier experiments and the second from later experiments, 
However, thase relationships are average and do not apply in some cases or 
throughout the complete possible range of permeabi lity and poresity, tere 
again grain sise, with ite attendant factors, directly affects these 
relationships. 

In waking deteruinations of the values of permesbility, effective 
porosity and crein sice diaseter, the first two are now very easily and 
effectively carried out by means of siaplified, standardized methods whieh 
give fairly asourete results, in the case of average grein size diaseters, 
this detercination becomes complicated by the particle itself, in that the 
particle is usually irregular in aise end shape, and no known sethods are 
avelisble for defining such o particle in geometric torma, Thie probles 
becomes sore complicated and forsidable in the case ef an azgregete of ir~ 
regular particles, but hes been successfully solved by the ase of statisti - 
cal aethods for evaluating surface ares and voluse of the aggregates, 


Tieforences are Listed in the Bibliography. 
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Such zothois are based an the avasurement of the so-celled average grein 
ameter by relating the irregular partiele to an equivalent geometric 
figure such as a circle or a sphere from which a diameter can be calealated, 
Thase sethods detersine surface area ani volume from the number and 

weight of the particles, using se @ besia the theory of ¢imensiens which 
requires that the surface area be proportional to the square, and the 
voluse be proportional te the cube of the dimensions used. 

Yor the deteraination of particle sise to have physical signi fi- 
cance, that is, to relate particle size to jefinite physical propertios 
euch es surface or volume, the erithaetic or ,eometric mean and the median 
dieseter determinations are superseded by statiatioal diamsters. Several 
formulas have been developed and used? to designate these statistical 
dianetere in terms of various physieal properties such as sean dianeter é 
(Aength), mean volumes diameter, mean surface diameter, moan velune-surface 
diameter, and weight sean diameter, ‘the average grain diameters used in 
this study are sean volume diameters, baged on the average volume and de~ 
fined to be that diaseter whose volume divided inte the total volume would 


give the number of particles, If d is the average particle diameter and 
¥ ie the totel volume, @ relationship existe as 
Vad nt? 


where n is the musher of particles of a given site, tf 4, represents the 
volume of an average particle, and Sn is the total nasber of particles 
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il. APSTORICAL BACKOROUND 


Investigation of the flow of fluide through porous aaterial at low 
pressures woe initiated by “arey? in 1556, ‘These studies consisted of a 
series of experiments on the flow of water through filter bede from which 
carey obtained an empirical formule showing the rate of flow to be very 
nearly proportional te the pressure drop per unit lencth of porous aediun, 
In ite elesentary fora, this formula, 

v a ~(5/.) OPAL) 
imilcates that the fluid permeshility, K, of a porous material propertion- 
ately relates the velocity, v, of flow of o fluid of u viscosity to the 
pressure differential, 4), causing this flow over a length of OL. ‘Since 
QO? is measured between the outflow and inflow ends, it is negetive as in- 
dicated in the formule above. | 
| The proportionality constant K is considered te be a specific 
property of the porous material, espiricaily independent of the dimensions 
of the material, the pressure differential exerted on the fluid flowing, 
and the viscosity of the fluid. For this reason, Kk may be expressed in 
teras of other measurable physical properties of the porous aaterial, ouch 
as porosity and grain aise. Since the formulation of Derey's equation and 
from it, his law, sany investigaters and researchers have attempted te show 
relationships of porosity, permeability, and grain else. : 

In 1880, Seslheim’? first introdueed grain size inte the relation- 
ship by showing that the rate of flow was propertional to the square of 
the average grein diameter, ‘this would indicate a decrease in permeability 
as sani graine decrease in Sizee 
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the next investigator to utilise grain size in flow relationship 
formulas was Hasen* who, in 1592 and 1993, proposed the forma based on 
stadies of sand fliters, 

qe Kp (4g) (h/t) 

where K, is @ permeability constant, h/L iv the pressure drop through the 
filter per unit~iepth, and d, ie the "effective" diameter of sand particles 
havin; uniformity coefficients of lese than 5. “fffective® grain sise is 
defined ag the opening which will just pass ten per cent of the particles 
(hy weight) and uniformity coefficient is defined as the ratio of the size 
opening which will pass 60 per cent of @ sanple being screened, to the 
size which wlll just pase 10 per cont. hen the uniforsity coefficient is 
low, sa 5 oF less mentioned above, the particles are sore or less uniform 
in size, but as tho uniforsity coefficient increases to higher values, the 
particle sizes become widely distributed. 

in using Hasen's equation, if q ie expressed in cubic fret of 
water per day per squere foot of filter urea, dy in willimeters, the lose 
of head, hy in feet of water, and Lin feet, then the values of K, range 
from 1300 te above 4000 depending on the cleanliness of the sami, with 
the uoual iisits of K, for ordinary sunds ranging between 2300 and 3300. 
Hasen's equation has been widely accepted and utilised by sanitery ongincers 
and is of considerable historical interest since it represents the firet 
attespt te recegnise the isportance of partiele sise and a sethed for its 
satisfactory representation, avia and Wilsey? retested Hasen's formule 
with the following results: 


as 2300 4,2 (SFB) nj) 


t # 10) 
ani 4 2 9080(4y,)°(6/u0)° = melas 
where T ie the teaperature of the water in degrees Fahrenheit, dy, is the 
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screen size in alliiasters through which just 3 per cent of the particles 
(oy weight) pase, and f io the porosity inisoduced inte the second formals. 
Glichter,” in i999, published the accounts of his theoretical 

analysis of fluid fiow through an ideal homegeneous packing conposed of 
epheres of unifersz sise in which he was the first to intredace the effeet 
of packing as @ factor influencing permeability, which cave rise to ex- 
pressing the average pore area in terss of the diameter of the sphere, 
Slichter derived a modified Poiseullle relationship in his equation 

@ = 10,2hd*b/K ul 
where q is expressed in oc of water per secon!, A ie the crose<sectional 
area of the peeking in square centimeters, h is the pressure differential 
in centimeters of water, 4 is the dismeter of the spheres in centimeters, 
K, 19 & packing constant dependant on porosity, 1 is the viscosity of the 
fluid flowing in poises, ani L is the thickness alony the direction of 
flow in centimeters. in ters of Sarey's permeability, Slichter's formala 
would read 

K = 10.24°/a, 
For various packings, tebulatiens have been made of the different values 
of the packing constant k, shown og @ function of porosity for porous 
material packed with spheres of equal diaueter, 

Wavis and Wilsey” developed a formula 

1/t, = 0.05(¢/40)"9 
for use in the deteraination of the values of K, where f is the porosity. 
King’ recognised the importance of Slichter's forma by using it to eal~ 
culate the average diameter of irregular particles, It is impertent to 
note here that the average dlanster of irregular particles determined in 
this senner is the diameter of the particles such that if a1) of them were 
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of thie dleweter, the packing woul heve the same permeability as it 
actually hes at q given teaperatare and porosity. 

For the ultimate in permesbility, porosity, and qrain aise rele- 
tionships, it resained for Keseny,° in 1927, to develop and publish his 
iapertent empirical formule 

v «(Pe oP /a8"L)(i/t,) 
where C is a form facter dependent on the shape of the pores, # is porosity, 
g io the acceleration constant, O? is the pressure differential in gram 
per square centineter, u is the viscosity of the fluid flowing in poises, 
6 te the apecifie surface aven in square contiantere per cube contineter, 
& is the thickness of the porous sedium in centiaeters, L, is the Longer 
path of flow in centimeters through the thickness L, and L/L, is a redue- 
tion factor for the relationship 4 P/L, Since § » 5,(1 - 9) the shove 
equation ean be rewritten in teras of 5,, which is the specific surface 
area por unit volume of satvix aaterial, on™, as 

¥ 2 (Pg ar/as,*(2 « #)*L)(L/i,) 

Unaware that Koseny had developed the aforementioned relationshlg, 
Pair and Hatehs? in 1933, developed an alucatidenticel formla by follew- 
ing @ line of reasoning siailar te that pureved by Koseny. This equation 
M | 

vs Par/sala « pis,” 
Tt ean be seen that the factors (C+L/i_,)> in the Koseny equation have 
been replaced by the constant $, believed by Fair and Sateh to apply to 
normal, rather compact, unconsolidated perows asterial, 

Carman”9et) 15 1937, presented an improved derivation of the 
Loweny and Fair end Neteh equations, Using as @ basic the general equa- 
tion for flow through pipes 
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v3 gn® A P/eal, 
whe re % ie the average velocity, ¢ is the aeceleration constant, m is the 
hydrwalic radius, defined as the volume of fluid in the pipe per surface 
area in contact with the fluid, 4? is the pressure differential, k, is a 
shape factor ani is the reciprocal of Keveny's form factor ©, u is the 
visoosity of the fluid flowing, and L, is the length of flow path, Carman 
set about to substitute or modify the equation feetors so thet the equa- 
tion would represent flow of fluid through porous media. Yor v he sube 
stituted v/f, later to modify it by the relationship L,/. for flow over 
the flew path, Ly greater than the i dimension of the medium, Then he 
substituted for a the relationship of #/ig(1 - #), end applying the redue~ 
ing fector i./i, to the expression 4 P/i, he arrived at the following 
rearrangement: 

veh 02/052 » #)"a)lu/iy)” 

in a series of experiments using a very wide range of particle 

shapes in unconsolidated media, Carsan found that the paraseter (0+L/i,,)"* 
of Koseny's equation and the paranoter k,(i./i.)” of his equation both 
averaged sbout 5, in the modified Loseny-Carman equation these parameters 
ere expressed as k, with an assigned value of 5, shown by 

vs 4 r/as.te(a ~ )*L 

ayinie,”” in nis paper on the historical development of the 

Kovseny<Carman equation, Jiseasees the more recent deterainations by ine 
vestigators of valuas greater than 5 for the Koseny-Carsan constant k 
when the formla hes been used in connection with measurements of flow 
through consolidated madia, Gince the k, shape feotor of the Keveny-larnan 
equation has been held by sost investigators to be constent, the higher 
values of k have been attribated to the (1y/i)* fustor which Sese end 
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grace) have called the tortuosity factor, T. Several deterninetions 

were made by jose end bruce for the value of k for various consolidated 
media of low permeability where values wore obtained as high as 100, 

Others coneur with ose and Bruges, all indicating that the constant for 
consolidated media varies to higher values and that these higher values 
are attributable to the tortuosity function. 

There has been some disagreement with the concept of increasing 
values of the koseny-Carean constant in connection with censolidated media, 
ile Ee Rose** velioves that errors arising from the use of 5 as a value for 
the constant will be far loss than those arising from ancertainties in 
other varishles, especially the porosity funetion, expressed as $°/(1 - #)* 
in the Kozeny~Carman equation, fie has shown, by plotting relative 
resistance versus per cent porosity, that at porosities Lower than 40 per 
cent there are increasingly different values of relative resistance to | 
flow for various equations which ali contain a porosity function of some 
sort, Sallevalle”” has shown a similar plot of reletive persesbility versus 
per cent porosity for several different investigators' formulas, all con- 
taining @ porosity function. ‘This also indicated that st porosities lover 
than hO per cent the porosity function docs not resain constant. It has 
not been possible to test the porosity function /(1 - ¢)* in consolidated 
media because it is not readily possible te vary porosity ond keep the 
specific eurface area constant. | 

A recent moti fication has been made in the Koseny equation by 
Wylde and Spangler who have casbined the Koneny equation with properties 
of the capillary pressure desaturation carve to obtain the formula | 

K w (y*/2.57*9) bis as /e," 


where K is the permeability coefficient, y is the interfeeial tension, F 
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is the formation feetor, # is the porosity fraction, 5, is the fractional 
wetting phase saturation, and P, 34% the caplilary preseure, ‘thie equation 
introduces new physical factors which are more easily sneagureble in 
laboratory ietersinetions and consequently should further the continuing 
search for mere accurate and useful relationship formulations for fluid 
flow through porgus sadia. 
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iii, PROCEDURE 


Por the purpose of carrying out this study, date of experimental 
and calculated vaiues for permeability, together with aecompanying values 
for poresity and average grain sise, wore selected from Fettke's*?* 
atudies of the Gradford (11 Pield, lredferd, Penneylvania, These dete 
hed been derived fros test determinations made on representative core 
sanplea obtained from various sections of the field at random depths in 
the Sradferd Third sand. ‘the iredford Third sand is an excellent repre- 
gentative of a congolidated porous sedium whose ceneral characteristics 
are low permeability and porosity and quite fine rain aise with irregular 
ohape ond aise distributien and a mxierate to high amount of 1i thi fleation, 
Sinee preetically ali cil and gas reservoirs are located in conselidated 
rook formations, a representative sanple of this type of porous aediun 
was dosen for application to this study. 

49 @ first approsch te the problem, an investigation was made of 
permeability, porosity and grain else relationships when applied te an 
equation modified by the author fron one developed by Trexler and Baus”? 
from the Poiseullie equation of flew through pipes. Thies equation was ex- 
pressed by Traxler and Baun ag: 


Bg 3 (32qai/aga ?)* (2) 


where D, ia the effective pore diameter in centimeters, q is the rate of 
flow in qubde centimeters per second, 1 is the viscosity of the fluid 
fiowing in poises, lL is the length of flow path in centimeters, A is the 
crose~ssctional area in square centimeters, # is the porosity, and 4? is 
the pressure differential in dynes per square centiaster, It can be 
readily seen that the function ql/AOP of this equation will be equivalent 
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te the peraeability coefficient Kk, in darcies, when the unite are changed 
to those of the Saray equation, ‘Thus, K « (qi00s./A9.8692 KX 107 OF) or 
945692 2.10 & 2 qul/An?, Sadstitating in equation (1) above gives: 


a, S629 X wp? {2) 


where K is the permeability coefficient in darsies. 

Table I, Appendix I, has been presented to shew calculations which 
resulted from the application of the permeability and porosity dete selected 
for this study to equation (2) abeve, solving for D,, the effective 
average pore diaaeter, Kutting’? nas stated that the effective average 
pore diameter is roughly one-fifth of the average grein diameter, Figure 1, 
Appendix If, shows « plot sade from the date in Table I of d, the average 
grain dlemster, versus 0, the effective average pore diameter, resulting — 
in @ straight line whose slope, which represents the ratio of d/%,, is 
approxinately 25. ‘The value of this ratio is coneiderebly greater than 
the one of $, (3/S4 = 0,), proposed by Butting., It should be explained 
here that only the results of tie caloulations using Pettze's 193, data 
were used in plotting Figure 1. ‘The resulte of the calculations using 
Pettke's 1933 date prodused no conclusive plot. 

Ag & second approach to the problem, an investigation was made of 
permeability, porosity and crein sise relationships when applied te a 
formulation derived by equating the KozenyCarman and the Darey equations, — 
This second study wag made on the basis of the specific surface area, 

Sq, of porous medium per unit volume of matrix material, Sy taking the 
bagie Koseny-Carean equation solved for 9," one obtains: 


2 Oe (3) 
*o ® a = 8)" qu 
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it can be seen thet is equation (3) on the preceding page the function 
(AA P/qud)™ i9 equivalent te the K of 2arey’s equstion when the unite of 
and O ? are changed from polees to centipolecs an! from grams per 
square centimeter to ataespheres respectively. ‘Thus, K = ql0GuL/ 

ris (2b P/qst), Substituting in 
equation (3) gives: 


(A9.678 x 10 VAP) or 9.670 x 10 
2 (&) 
& — 
o * Fa - 989.678 & 


Since the seceleration constant, gy is 930 centineters per second squered, 
equation (4) wbeve can now be resolved as 


Sy 2 hsare fs oi? ‘s) 


whieh in the equation used in this second approesh to the problen water 
etudy, | 

Table Il, Appendix I, hes been prepared to show caloulations re- 
eulting from the application to equation (5) above of the permesbi lity 
end perenity data selested for thie stay, 1% hes been chow ty Serden”” 
that S, = 2,/éA for irregular particles where Z,/ly i9 the ratio of the 
eurfece to volume factors ond d is the average grain diameter, In this 
ratio, 2, 2 5/d,*, where 5 is the surface ares and 4, is the statistical 
particle diaaeter based on surface arsa, and Z, = V/dy, where Vis the 
volune and d. is the statistical particle diameter based on volume. ‘The 
ratio 1,/i, has been imifeated by Pair and Hateh? to be a useful measure 
of particle shape ani has been given @ value of 6,0 for spheres, 6.1 for 
rounded particles, 6.4 for worn particles, 7.0 for sharp particles, and 
7e? for angular perticles, Pigure 2, Appendix II, shows a plot made of 
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the values in Table 11 of &, versus 1/4 which roughly fors straight Lines 
with a slope of 143 for Line 1 using Pettke's 193) deta and a slope of S$? 
for line 2 using Fettke’s 1935 date, The slope normally should represent 
the values of the ratio 2/2, which iie somewhere between 6,0 and 7.7. 

As @ third approach to the problen, an investigation was made ef 
permeability, pewosity and grain sise relationships when applied to 8 
formation derived by equating the foseny-Carman ani the “arcy equations, 
as previously aentioned, except that the formation was modified for the 
use of average grain size in lieu of specific surface area, Since , is 
proportional to 1/dy 8,” can be replaced by 1/4” in egastion (ii) abowes 
There is also a new constant, k', introduced inte the equation whieh is 
equivalent to the function 1/(2,/z,)°5 where 5 is the Keveny=Carnan eon 
stant and 2,/i, is the ratio of the surface to volume facters, ‘The re~ 
Salting equetion is expressed as: 


3 
(2 = $)*9.678 x 
By oubstiteting 980 centimeters per second squared for g in equation (4) 
end rearranging to solve for k*, the equeticn can now be expressed ag: 
k's k= 9)* (7) 
1.013 x 
Table III, Appendix I, hae been prepared fron caleuletions derived 
by applying the date used in this study to the factors of equation (7); 
namely, K/d? and 2.023 x 10° /(a ~ 9)", the values of these factors were 
plotted as /é? versus 1.023 x 20°99/(a ~ ¢)*, show by Migure 3, 
Appendix Ii, to be a straight line with a slope of 0.000L052 in the case 
of line 1 using Fettke'’s 1934 data and a curve with « variable slope of 


as 


eons AR dowe emot edyuor moite SL anerer ag? © To adder ot conten one 

QE So 'dyale a tne sont MEE eteateot yasns £ ont det GAL to eqote a Able 

seavesgyt Mbuoite Yimirina eqnde edi sagas URL etodioat poten § ati oo 

Be Galena este tay AR eds ote Ye sensor eas 
ae ON ne hinting at os Contras Dada © at 


ia'y biaty Cm 

me Ca a Fas A, 

e y > may 4, nf . ao > pe ; ~ be ad er ‘ist , faba Lay 
Ys Lx “ese Jia 


: arate hy bene® - mpanie ks sa tring ‘we oa 2 ee RAS 


na we es 


= met 


deste 3s Hotan dome Ga aehilics 4 42 0 sil aliens 
To eqole eidsiray © idle wrwe @ bow osb MCS eadddet gaten £ ents te 


unknown value in the cage of line 2 using Pottke's 1936 data. Line 2 
could not be expressed as @ straight line by repletting on semi-log or 
leg-leg graph paper. the slope of line 1 of Figure 3 ic equal to the 
constant k* of equation (7) on the preceding pages 

As previously stated, k' is equivalent to the funetion 1/(z,/2,)*5, 
Uaing the value of the slope determined for line 1 of Figure 3 for k*, the 
value of 2,/i, as obteined fron the relation Zg/ty = (1/k'5)? 48 35. 
formally, values for i,/Z lie between 6.0 and 7.7 a8 previously aontioned. 

Im the two preceding investigetions, the porosity function, 
#/0. ~ g)*, sathough At varies with the value of porosity, is considered 
& constant factor in that it reteins the same fiwed relationship regard- 
less of the source of data, For s fourth appreseh to the problem, an in 
vestigation was made of the permesbility, porosity and grain sise relation- — 
ships when applied to a formation developed to expreae the porosity func- 
tion a9 @ variable, dependent upen the sourve of data, with the persee- 
bility and grain size being constant factors regardless of the source of 
data, ‘Thus, the porosity function would be 6 constant factor only for « 
particular reservoir formation or even a particular section of the 
reservoir formation. 

To develop this fonmlation, the basic Koseny-Carnan equation was 
recast to include particle diameter, giving: 

‘ ra 

whese units have been previcusly defined, Solving for the porosity function 
gives: 
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Sinee the function qiL/AD?, when expressed in terms of the “erey equation, 
is equal to the permeability coefficient i of Carey's equation, eubetita- 
tion of 9.678 % 107K for qul/ADP can be made in equetion (9) above civing 


GQ =) n g k* 


Aseuming that the constant k'! would be reflected in @ recasting of the 
porosity function, @ new expression of this function can be substituted 
for kP/(1 - g)* in equation (10) above giving 


(#)" e & / 1.003 x 108s® (aa) 
or, by rearrangesont, 
| K = 1013 x 20° a (g)” (12) 
where n is a constant applicable only te 4 pertieular formation or section 
of a fornation, or & porous medium, 

Table IV, Appendix I, has been prepared from caloulations derived 
from the application of Fettke's 193 data to the functions ef equation 
(12) above. Figure 4, Appendix Ii, is a leg-lez plot of these functions, 
«/42/.,013 x 10° versus g, which forme @ straight line with « slopes of 6, 


oy a = 205 t/ah aca 1 0h 


then n ¢ 6, the slope of the straight line. 
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iV. GOMSENT AND ULSOUSSION oF BZSULTS 


Before discussing the results, it should be pointed out here thet 
the values of porosity data used in this study represent absolute 
porosity. However, Fettke® has shown that in a eajority of eases the 
difference between ebaolate porosity and a Lower value for effective 
porssity in the Draiferi Third sand is but o matter of 6,1 ser cont 
averages. There are & fow isclated cases of @ deviation slightly greater 
than OG.) per cent up to about 1.0 per cent. For purposes of this stady, 
the values ef effective porosity sere assumed te be the same as those for 
absolute porosity, This fact can explain some of the siight deviations 
of the pointe used to plot the various lines ani curves shown by Figures 
1, 2, 3 ond &. Other devietions of the points can test be expli<ined by 
the facet thet at the time the data were generated, the securasy of the 
éetersinations of values for permeability, porosity and grein size was not 
nearly es good as that of the present day, due to greetiy improved and 
Standardised sethods for deteraining such values. 

the results of the author's firet investigation clearly show that 
the average pore diameter in a rather Aine=grained, consolidated median 
decreases considerably from the average for unconsolidated media, stated 
by Batting”? to be roughly one-fifth of the average grain diaseter. This 
ean be best explained by the non-uniformity of particle size and shape, 
the fine grain sie, and by the degree of 1ithifiestion or conenting 
found in consolidated media, It is considered by the author that the 
latter exerts the createct effect upon average pore dlemeter since cerent- — 
ing to any extent tends te restrict, reduce and plug the pore channels of 
the rock, To utilise this mathod of correlation of peresability, porosity 
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ani grain eise, it would be necessary to ieotermine the relationship be- 
tween averaje affective pere dlaseter and average grain diameter for each 
particular reservoir or section of the reservoir, 

im the second approach to the prebles, the results of the in- 
vestijation indicate that the vaiue of the Loseny~-Carman constant for 
consolidated sedia is auch hicher than the assigned value of 5. fy way 
of exasple, preaume the consolidated sedium to be composed enerally of 
sharp particles where the valae of the ratio of surfsce to volume factors, 
Sg/ty) ia 760, 1% haw beon previously shown that 5, @ 2,/i,4, Using as 
@ selevtion from Fetthe's 193) data, 4d = .0093 centimeter, i » .0034 lerey, 
and @ = 2145, it is possible, by resolution of equation (4) of section IIT 
of thie study, to soive for a value of the Leseny-Carnan constant in lieu 
of the value of 5 normally used. 5, is resolved as 7/,0093 or 753. 
Solving for the unknown doseny~Carnan constant using the above values in 
the equation gives a value of approxinately 220, by using som of the other 
sets of values for permeability, porosity, and grain sise from the 1934 
date to solve for the foxeny~Ceraan constant in the save sanner as described 
above gives valuep for the constant in the senerel vicinity of 226 with 
some variations, the results of similar application of Fettke's 1933 data 
were not as conclusive, with the values showing a meh wider variation, 

3 to 107, for the Kozeny-Carnan constant, 

4a in the second approach, the third approach to the preblem alse 
indleated that the Koseny-Caraan constant for consolidated media is each 
higher then 5, it has been previously stated in Section Ill ef this 
stady that the elope of Line 1 of Figure 3, Appendix IT, was 0,00010S2 and 
equivalent to the function 1/(Z,/ly)"k where k 49 the Keseny~Carnan 
constant. if, for example, it is again assumed that i,/Z, hee a value of 
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7eO for this medium, it is peasible te aclve for a velue of other than $ 
fer the Kexeny-Caraan constant, in this cuse, the value was resolved as 
19% for the congiant. ince @ small deviation in the drawing of line 1 
of either Aigure 2 er 3 results in 4 change in the value of the alope of 
the ling, it can be presumed that this feet osused the di fferenee in 
values for the (osenyeCarnan constant determined by the two investigations, 
Therefore, 1% can also be presumed that in the cave of the 1934 date used 
the averace value for the igneny~-Carean constant lies scnewhere between 
194 and 220, 4 sialler suceese was achleved in using Yettke's 1939 date, 
it ig believed that the 1939 data were not traly representative of the 
reservoir sand oer were they in theaselves complete enough to peruit dras- 
inj ary coneiusions. | 

the fourth epproagh to the probles was an attespt te show the 
porosity function as a variable rather than @ constant faehor in liga of 
@ variable for the Loveny~Carman constant as shown by the secon? and third 
investigations, 4n application of the 193) data to equetion (12) of part 
TUE imileates that the » of the porosity function, #, is approxiaately 
eqaal to 4, ae deterained fron the slope of the line plotted in Pigure hy 
Appendix il, with higher and lewer variations in sou cases, Since the 
porosity function is an arbitrary selection by the author, it follows 
that thie particular function, being an exponential, could have a variety 
af foras all of whieh would, when solved numerically, cive substantially 
the sane values, @ was adopted because of ite singlicity if for no 
other reagon, Although this aethod shows thet the peresity function can 
well be a varlable factor, it is believed that the worksof Sese and Bruce? 
and others have contlusively shown that the porosity function resains 
constant wile the oseny-Caraan constant becomes a variable factor in 
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eonaolidatet aedie, %e attenpt was sade in the fourth investigation te 
make use of Fettke's 1938 data because of provicus inconsequential 
results 

im ali of the author's four investigations, it has been show 
that the procedures, equations, ani constants applicable te the empirical 
formaletions derived from the studies of unconsolidated percus media ds 
not hold true in the case of consolidated perous aodla, However, it siso 
hae been shown that for a ,iven or particular consolidated porous rock 
formation, the eapirical formations apply when the norsal constants are 
changed to those values which fit the sitastion, ag im the case of this 
stady of the consolidated Sradferd Third send. ‘The evidence (oes not 
conclusively prove whether the porosity function should be « variable or 
& constant factor in persesbility, peresity and grain size relationship 
formulations. ‘uch mere invrestigution ani research wust be carried out 
alien; these lines before @ conclusion can be reached ag to which theory 
is correct, although ao mentioned in the preceding paragraph present 
eviience pointe to the porosity function being a constant factor, It is 
very possible that beth the porosity fanction ani the loseny-Carsan con~ 
stant aight be found te be variables in the ease of consolidated porous 
muectia, 
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The fleld of email particle statistics and asagurevents is exten« 
sive, comlex and relatively unknown, Yor these reasons, the stady of 
particle size in relation to permeability ani peresity of porous media is, 
of iteelf, extensive, complex ani relatively anknown, especially in the 
ease of consolidated porcas media. the author agress with tuskat@ who 
stated, in effect, that the use of snall porticle atatietics ani measure- 
nents in ieterninations of permeability was so complex and unresolved that 
che cresently utilised methoie of detersinin: permeability and porosity 
fer outweigh any aethod using grain eine or ite atteniant properties to 
detersine same. However, the suthor feels that it ie eusentially im 
portent to intensively stady the field of particle sise in order te qein | 
a thorough knowlede end anierstanding of macroseepic fluid flew in perous 
media in the light of grain aise and its attendant fsetors of pecking, 
unifornity of sive and shape, and degrees of lithification, all of which 
affect fluid flew in such media, 

Through an intensive study of thie nature, 1t is possible that 
some day a forwalation will be developed which will truly hold uniter a1) 
condi tions end in any circumstances whereby the physical proverties of 
permeability, porosity, and crain sise with ite attemiant faetors, can be 
successfully correlated. ‘the four methods of correlation eapleyed in 
thie etady can be used in a restricted degree when applied to the cen- 
seliiated porous sedia of particular reserveirs or of sections of the 
reservoirs only efter considerable researeh to detersine the proper 
constante applicable in eech case, in the same vein, and of considerable 
interest, is the method developed by liyter”* for makin; pernssbility 
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anenaremate without the use of core sanples by weans of an eopirically 
developed formula, good only for one particular sani, which uses only 
jeterained weluhia of sieved sani particles which have been passed 
through a series of U. 5. Standard eleves, 

_ Another possible investigation or sethod, suggested to;, but un- 
tried by the author, would be to plot the results of sieving data fron 
@ core sample (per cent by weight of sanple through and on each sereen 
sise) versus the average soreen aise for that per cent by weight. Pron 
the resultant carve, select a "representative" per cent value anid deter- 
aine the ¢gorreapomding average sereen sive or presumed averaze grain 
diameter, Using equation (2) of Section [11 of this etady, deternine 
whieh per eent by weight value and corresponding sersen sixe or presamed 
averace grein cise can be used to satisfy the equation, ani then deterszine— 
whether or not this sane per cent value can be used in connection with 
@ll core sacple data from cores taken from the sane formation or section 
of the formation. 

A #amaation ef the rewalte of this study does not reflect the 
original purpose of the stady, I4 wos thought thet on evaluation could 
be sade of the interrelationshipe of the physical propertics of porous 
rock, permeability, porosity and grain sise, as well og @ modification of 
existing formulations used te express these interrelationshios of physica: 
properties, Ry using date from bichiy consolidated porous sedia, it was 
hoped that if a modification were developed for effectively «expressin; the 
interrelationship of persesbility, porosity, and grain sise, it would hold 
quler ang chveunstenses ent conditions fev ony type wethun, oi ther cane | 
Solidated or anconsolidate’, owever, the complexities of grain sise 
stady, the dearth of sufficient and accurate data, ani the apparent non- 
eonforvity of apolieation of medified formlations to consolidated porous 
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media data sade such endeavors impractical an’ impossible, 

It ie believed that the regulis of this study, although meager, 
can contribute some benefits in a better unlerstanding of the problex, 
the ratio of pore volume to jrain sise diameter has been shown to decrease 
markedly with an increase in cementation and a correspending decrease in 
grain size euch as found in @ consolidated perous media like the Bradford 
Third sand, ‘the value of the Aozeny~Carnan constant has been shown to 
inerease to very high values in consolidated porous media, This is preb- 
ebly due to the large increase in values of the tortuosity factor which 
becomes greater as the rock beeemes sore consolidated because of an ine 
crease in the effective length of flow path. frobably the fost iaportant 
contribution of this stady io the realization of the importanee of the 
Koseny-Carman equction and ite aodtifications when applied to problemas of 
fluid flow through porous media, It holds great possibilities, when sub- 
ject to the proper modi Meation, in ice applicetion to stadies of this 
type ami gives promiee of eventual sclution of the sany probhems of fluid 
flow still faging the petroleum industry. 
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TABLE 7 
Fettke's 1934 “ate (kK ie Experimental) 


(t) 


Grain Sise Permeability Porosity 4/# (8/8) constant Xe 3 aa/iy! 
r & (Pargies r cent | x 14 (om) x 109 


afX92 20005 16.8 «OGh50 .0678 5,619 3.81 
0092 20005 12.5 +0Gk00 0633 * 3256 
20093 0ORy Lied 202340 .1530 ° 3.60 
009% 2082 15.3 205290 42278 * 12,80 
20699 0012 130% 00895 20946 ° 5.32 
20106 e009 15.8 082326 «.1523 ° 8.57 
OL09 15.2 eO2LKO 1562 « §.78 
Pry ees) 20039 16.2 02410 41552 * 3.7% 
Pays Ts) OL 39 202 «06590 .2625 * 7S 
eOL2L 0213 21.6 eO987G «5 R12 « 17.00 
OL23 +0366 22.7 e802 * 22,50 
e132 00323 21.7 23860 " 
«0226 Lie? eOL?7O 41330 a Tek? 
+0139 -OOLT 15.0 ° «1063 * 5.97 
Fettke's 1934 Sata (Ek is Caleelsted) 
20092 «000k oan 20608 2=S ss § 49 3ek2 
e092 20080 12.65 205% bd §_02 
pome . Lied 0660 41288 . * 7023 
15.8 s02h00 21549 ° be Th 
20099 OOS L304 +0200 1095S ° 615 
20106 20085 16.8 202680 .1637 e 9-280 
«09 20031 15.2 eO20h0 41425 ® 5.02 
O26 20044 202720 .16N9 a 9e26 
GAO oO Ss 2002 207630 2762 * 15.50 
0121 20225 21.6 = 20800 +3285 " 
0028 20293 2207 013130 * 
eOL32 20232 2he7 030600 43256 * 15.30 
2135 20026 Lie? eQL770 41330 ° Tolee 
20139 90029 15.0 oOL930 «41389 ° 7280 
Pettke's 1930 “sate (EK is izperimental) 
20041 000278 #02090 lhbhG 5619 5.12 
900R2 0000K7 12,1 200389 0624 . 350 
e0OU2 /OOL31 12.6 $.O1M10 .i1187 ad 68 
20042 «00550 lied 003840 51960 * 11,00 
2002 s0OV12 205930 a ® 13.72 
200b4 «00534 ihe3 20370. ® 10,36 
20047 eO2752 1507 eAT530 44137 ° 23.50 
20051 02775 15.7 217630 44199 a 23,60 
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TABLE Ti 
Fettke's 1934 cata (kK ie Secerisental) 


ONDE lak ee ee 


rT} 3) 4 9} rss TS} ; 
Grain Sise Permeability ferocity L+-¢ (buy Seegi fle 1/a 


a fom) K (Sareies) / (per cont) furfage 3, 
( curl} . 
20092 26005 16,8 0392 ike 70 5000 103.6 
0092 29005 12.5 e875 «61S HE — 108.6 
e093 003% Lied 0855 e533 LOT 
00% 22082 15.8 2342 in 39 3700 106.3 
eOOP9 eOG12 L3eu 0866 7350 201.0 
eOL06 GOS? 16.8 «532 %e57 5970 Piao ds 
2009 0039 1562 og Geil S166 8 
Ohl 0039 16.2 538 SekS S6L0 Fl. 
eOLlO «0139 20.2 a 345% 4340 PLO 
eOh21 0213 Z1.6 2 78% 3045 3950 $2.7 
oL28 > 2267 o?73 2 eu? 3290 73.2 
e032 «0323 Phe? 07830 2459 3230 7508 
O35 . Lite? Ve52 5525 Mek 
«UL39 SOL? 15.0 350 9.39 ThGo Tie? 
Pottibe's 1934 cate (K ie Caleulated} 
Cove 20004 10,8 e392 «= 16,43 8960 108.6 
e010 12.5 e575) 3=61318 TLS 103.6 
0093 20024 lied 2355 Ve?? $930 107.is 
0% 20533 25.8 odu2 5456 106.3 
0099 20015 136k — @fS6 Veled 6580 1aa.o 
00166 oOOkS 16,8 +832 G12 5630 Piso . 
2009 e0OR2 152 2848 7.00 5650 71.8 
Pes) Sts) 16.2 833 527s Vid 
LLG »SL5h 2002 o 738 4120 71.0 
e121 20225 216 78h = 3.0 sao 82,7 
0828 0298 2257 773 2.76 3685 781 
oh 32 «O31 21.7 «783 3330 7508 
e035 «0026 Las? e853 0=— (7452 6525 Thed 
«139 «0629 15.0 «350 Tel §700 Tie? 
Pevtke's 1933 Sata (K is Experimental) 
2CORL * 13-3 2867 Sei 770 Zidane 
«0042 -0O0h7 12,1 2379 9960 233.0 
oOGu2 . 12,8 «872 Selk §S66 236.0 
+0002 +Q0550 1423 ed 5.10 3880 235.0 
20Qk2 e0O9L2 152 ° &.cd 3300 238.0 
(Qld 200534 4a} e857 = SB 3390 227.0 
20047 202752 4567 2039 2000 212.5 
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Pettke's 1934 date (Kk is experimental) 


Grain Size Permeability Poresity ,. a-g* n/a 2/6)" 3 
20092 90005 10.8 800138 =. 795 5.971 176000 
00092 20005 12,5 0089S =. PGE 5.9 259000 
20093 20038 MaeS 200305 «= TR 39.30 423000 
200% «0062 15.6 0039, «=. 708 953.20 S64000 
O99 20012 e002 «=. THO 2 326000 
00106 20039 146.8 eOOh7, 2692 OF4OOO 
2009 20037 1592 #00380 720 493000 
2810 0039 eO0K25 «70h 32.20 613000 
20116 20439 2002 0082 «694634 «115.00 1920000 
20122 20213 21.6 e01007 = GS 1650000 
20126 20366 2247 01180 6.597 «6223.00 §46—2ono60 
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Fettke's 1994 data (K is caleulated) 
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Fettke's 1933 date (Kk is experimental) 
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K (Sareies) 4 (em) 1,013 x 10° ¢% (per cent) 
20005 0692 5.9% 5.83 x 1078 10.8 
20005 0092 SoA 5,83 12.5 
20034 20093 39930 38.70 UaeS 
20062 009% 92.60 15.8 
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«0139 2OLLO 115-00 113.50 2002 
20213 20421 146.00 iil 200 21.6 
00366 20128 223.00 220.00 22.7 
00323 20132 165.60 183.00 21.7 
22026 #0135 Linn 0 1.20 1ae7 
0017 20139 8,80 8,68 15.0 

(Kk 49 caloulated) 
200k .0092 be 72 466 x 10° 1068 
+0010 oOO92 11.30 11.64 12.5 
2002, wy 27.80 27k weed 
+0038 « 4300 42.40 
«0015 20099 15.30 15.10 136k 
20045 Prat ts | 40.20 3960 16,8 
0031 20109 26.10 25270 1502 
200kk -O110 35.90 1662 
oOL54 eOLL0 123.20 126.60 2002 
20225 e121 152.00 21.6 
00298 20123 182,00 179.60 2267 
00232 20132 132.380 431.00 2167 
0026 «0135 14.30 Lbekd Use? 
20029 20139 15.00 14.80 15.0 
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